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O : ABSTRACT 

The effective recombination coefficient, Oefr, is refined for optically thin cases. 
Radiation pressure in Lyman-a and Hell Lyman-a is calculated in terms of 
^ ' the number density and the mean free path. Pressure equilibrium between 

^ ■ Lyman-a clouds and an adiabatically expanding inter-galactic medium is as- 

O ■ sumed, n^Tj = ricTc- Numerical models of isothermal and adiabatic expanding 

^ , Lyman-a forest clouds are presented, including evolving UV sources (QSOs), 

^ ! with various turn-on times Zon = 20, 10, and 5, and with = ^ in a matter- 

\ dominated Priedmann-Robertson- Walker Universe. These models lead to the 

' conclusion that the radiation pressure and QSO turn-on time are significant in 

determining the range of physical size, D, and neutral hydrogen column density, 
^ ' N{HI), permitted for stable Lyman-a forest clouds. 

■ This manuscript was written in 1989 and never submitted for publication. 



1. Introduction 



A great deal of effort is being expended on the absorption features of QSOs, especially the 
IGM clouds that comprise the Lyman-a forest. Although the method of confinement for 
these clouds is still questionable, current research provides a detailed list of physical char 



acteri stics, including: physical size, column density and temperature, see ICarswell et al. 



(|l987l ). Sargent et. al (1980) and Chaffee et. al (1986). Typical values for the physical size 
of the clouds is 0.4 kpc < D < 400 kpc, Sargent (1988) and for neutral hydrogen column 
densities is N{HI) ~ 10^^ cm~'^. In addition to these physical characteristics, the Lyman-a 
clouds exhibit a distribution in redshift and neutral hydrogen column density of the form 

P{NHO,z)^N-^{l + zr , (1) 
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where (3 ^ 1.7 it 0.2 and 7 ~ 2.3 it 0.4, and the expression is vahd over the range of column 
densities lO^^ cm'^ < < 10^^ cm-^ Carswell (1988) and Sargent (1988). 

Numerical models of Lyman-a forest clouds showing that the upper limit in column density 
for the distribution is due to radiation pressure are presented here. We adopt the quasar 
evolution model proposed by Schmidt & Green (1983) with turn-on at z = 20, 10, and 
5, and follow the size and density evolution of the clouds that were stable to radiation 
pressure at the time the quasars formed. We show three models. A, B and C, by varying 
the magnitude of the background UV flux within accepted uncertainty. We also calculate 
an effective recombination coefficient, Ocfr, for determining the fractional ionization of the 
cloud in the case r << 1 as a prelude to the model simulations. 



2. The Effective Recombination CoefBcient 



The radiation field is normally separated into two parts, a "source" part, resulting from the 
background radiation field, and a "diffuse" part, resulting from the emission of the ionized 
gas, Osterbrock (1988). The ionization equation for a one element gas is then 

/•oo t(s) _|_ Ad) 

Nhi / " T' " (^udu = NeNpaA , (2) 

where a a is the total recombination coefficient. The "source" part is assumed to be a power 
law, Fi, ~ z/~". For a plane parallel cloud of thickness T^^, the "diffuse" part is given by 

where P = hi^Q/kTe, Mihalas (1978). We then solve the ionization equation to obtain the 
fractional ionization with optical depth. For small optical depths, tq < 1, the fractional 
ionization lies between case A and case B, and approaches case B for large optical depths 
To > 1, Figure m 

We desire a generalized numerical form of the efficiency factor that is valid for a variety 
of optical depths and electron temperatures. Therefore we solve the ionization equilibrium 
equation, including both source and diffuse terms, for different optical depths and temper- 
atures. A simple function of optical depth, optical thickness and temperature is then fit to 
the solutions as the efficiency factor e(r, T, Tg) 



The general ionization equation is 



t(s) /-oo j{d) 



I ' Iw^"^^^ / "hv^"^^] = '^e?^paA(7'e) • (5) 
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Fig. 1. — The fractional ionization, ^ = nni/riH, as a function of optical depth, tq, for 
a semi-infinite medium with = 10^ K. Case A uses the total recombination coefficient, 
Case B ignores recombinations to the ground state and Case E is the effective recombination 
coefficient. 



The ionization state of the gas does not depend on the nature of the source of the ionizing 
radiation, only the number of ultraviolet photons greater than threshold. Therefore we 
take the source of ionizing radiation as an active galactic nucleus (AGN) characterized by 
a power law dependence on frequency throughout the ultraviolet portion of the continuum. 
This power law is characterized by spectral index a and the ionization parameter F 

~ z^-" , (6) 

r ~ ^ , (7) 
riHC 

where Qo is the number of ionizing quanta per unit area per unit time and nn is the total 
hydrogen density, uh = np + uhi- 

The above form for the source term allows the first integral in Eq. (0) to be done analytically. 
Typical values for a, F and Qq are 3/2, 10~^ and 3 x 10^^ cm~^ s~^, respectively. Then 
the number of photoionizations per unit volume per unit time due to the "source" radiation 
field is 

S = UHI ^^-±la,e"^-^du , (8) 



nniCTo 



oo 
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where x = v/vq and Qq = Qx ^/'^dx. The substitution y = tqx ^ allows the integral to 
be evaluated 

n=0 

which gives the number of photoionizations per unit volume per unit time at optical depth 
To into the cloud. 

For To < 4, the sum rapidly converges and in this form is easily used. For tq > 4, ^ = 
0.40157rJ"^'^^^^ is an adequate approximation. The ionization with no "diffuse" part is 

°° r" 1 

ra///Qoo-o V!(-l)"-T7^ — -^=nenpaB , (10) 

n=0 

which is easily solved by letting rig = Up, and defining ^ = nni/n-H- These substitutions 
result in a quadratic for the fractional ionization ^ 

e-42+'-^±i-ir^-i-^]+i = o . (11) 

\ as n n! 2n + 3 / 

\ n=0 / 



The ionization equation is now solved including the "diffuse" field by integrating Eq. ([3]) 
over frequency 



tW poo „—(3x roo 

- / f dx 

lUQ Jl X Jq 

where x = u/uq and 



roo i-oo „—px roo 

nm ^a,dv = f{T,)nl dx E2{t,)[1 - i{T,)f dr, , (12) 

Jun Jl X Jq 



Using the definition for the exponential integral allows us to write the full ionization equation 



as 



(1-0^ = e^D-ir^^ (14) 

OA r. n\ Zn + 6 

71=0 



f{Te)Em 



OiA Jo 



poo 

/ S2(ro)[l-e(To)]'dTo . 
Jo 



This equation is then solved for as a function of optical depth, tq, shown as the calculated 
efficiency in Figure [2j 

These same results can be obtained by considering an "effective" recombination coefficient, 
= a_B + ^{t, r, Tg) X ai, and considering only the "source" radiation field. The "usual" 
approximation for the efficiency has been e = ^e~'^. While this approximation works quite 
well for thick clouds, it does not accurately describe the situation for thin clouds as in the 



- 5 - 



Efficiency Factor for an Optically Thin Media 
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Fig. 2. — The efficiency factor e(T, T, Tg) vs optical depth for media of optical thickness 
T = 1 and temperature Tg = 10^ K. 



case of Lyman-a forest clouds, Figure [2j We present here a new form for the efficiency 
factor e(r, T, Tg) dependent upon optical depth r, optical thickness T and temperature Tg 

e(T,r,Te) = /(/+ (1 -5)tanhy)((l -tanhx) + (1 -tanhx')) , (15) 

where 

y = -0.3 log T , 

X = 1.25(logr + (16) 
x' = 1.25(log(r - r) + /i) , 
te = Te/lO^ . 

The expressions for /, g, and h for are given in Table [H For Te > 10^ K, the efficiency 
factor becomes 

e = inf[l, (1 - 0.007 exp(-0.00338te - 0.573T))] . (17) 

We used the photoionization code CLOUDY and our new model for the efficiency to repro- 
duce the results from van Blerkom and Hummer (1967) as a check of the numerics. 



- 6 - 



3. Evolution of the Background Ultraviolet Flux 



The number density of UV photons per unit frequency, n(t,i'), will satisfy the continuity 
equation 

dn . d f du\ 



_ + V(,.v) + ^^^,.-J=S„W , (18) 

where v is the expansion velocity of the Universe and Sy{t) is the source function for UV 
photons. Assuming that UV photons are continuously supplied by quasars, then Sy{t) is 
proportional to the number of quasars. Following Ikeuchi &: Ostriker (1986), we adopt the 
quasar evolution model proposed by Schmidt & Green (1983) with go = = 9, z^n = 20, 
10, 5 and 2.5. Then the mean intensity at the Lyman limit is 

At:J^^{z) =A-KhvTc{l + zfS,,^{'d)F{z,Zon) , (19) 
where we have assumed Sy{ti) i^^^ and F(z, Zq^) is defined by 

with t{z) as the fractional look-back time to the present age. 

We define ^^^^(O) such that A7tJ„^{2.5) = Att x 10"^^ ergs cm"^ Hz~^ x /, Ostriker & 
Ikeuchi (1983). We vary the factor / in our models from 1 to 100 (models A, B and C). 
The ratio Jy^{z) / Jy^{2.h) is shown in Figure[3]for the evolutionary models described above. 
Figure H] shows that the redshift range 2.5 < z < 4.0, the Schmidt-Green QSO evolution 
model produces a background that is nearly constant and in agreement with Bajtlik, et. al 
(1988). 



We let the IGM and the Lyman-a clouds maintain pressure equilibrium and consider two 
cases for the expanding clouds: the clouds expand isothermally, or the clouds expand adi- 
abatically. If the IGM is expanding adiabatically and the temperature of the IGM is non- 
relativistic, then 

n,(Z)ri(z) =ni(0)(l+z)3ri(0)(l + z)2 = ni(0)ri(0)(l + z)5 . (21) 

Then for isothermally expanding Lyman-a clouds, Tc{z) = Tc{Q) and the density, ric and 
radius, i?c, for a fixed cloud mass scale as 

ne(z) = ne(0)(l + z)5 , (22) 



Factor 


Te < 25, 000 K 


25, 000 <Te< 10^ K 


/ 


0.2510 exp(0.001518te) 


(0.4346 + 0.0191 log re)/2 


9 


0.9814 exp(-0.002208te) 


1.167 + (0.0035 - 0.0144 log Tg) log 


h 


0.3739 exp(-0.006472ie) 


log(3.3429 - 0.2794 log Te) 



Table 1: Best fit for the factors /, g and h. 
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Redshift History of Ionizing Flux 
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Fig. 3. — The ratio of UV flux at redshift z to UV flux at z = 2.5 (at threshold) vs. redshift. 



and 



Rc{z) = Rc{0){l + z) 



(23) 



For Lyman-a clouds that expand adiabatically, the density scales as 

ndz) = n^iO)il + zf , 
and the radius of the cloud follows the expansion of the Universe 

Rc{z) = RMi^ + ■ 



(24) 



(25) 



4. Radiation Pressure in the Lyman-a Lines 



These primordial clouds are optically thin in the continuum, tq ~ 10~^. However, they 
are optically thick in the Lyman-a lines of hydrogen and helium, Tq ~ 10^, and therefore 
radiation pressure in these lines can be significant. The radiation pressure is 
where is the number density of the Lyman-a photons. The number density is found by 
balancing production and destruction of photons, Mathews (1976) and Elitzur & Ferland 
(1986). Lyman-a photons are produced from recombinations to excited states 



production = Ug^UpaB 



(26) 
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These photons are destroyed by repeated scatterings until the photon escapes from the 
cloud. The number of photons that escape per unit volume per unit time, nesc, can be 
found from the number of Lyman-a photons divided by the average time a photon spends 
within the cloud. This time scale is simply the mean free path divided by c. Therefore 

destruction = tIqc/L , (27) 

where L is the mean free path. Bonilha et. al (1979) provide an analytical expression for 
the optical mean free path, Lq. To convert Lq to mean free path, we must divide by the 
line absorption coefficient, = ^hi^o- Equating production and destruction and solving 
for Ua 

LoUeUpaB 

Ua = — , (28) 

and for the radiation pressure 

Hell Lyman-a is treated differently since A304 photons may also ionize H and therefore have 
an additional destruction mechanism. Bonilha et. al also give an expression for the optical 
mean free path in the presence of absorbers. Their R, which is the ratio of optical depth of 
absorbers, in this case T912 and the optical depth at line center, t^q^ is further reduced by 
the photoionization cross-section of hydrogen at 3 and 1 Rydbergs, a^/ai. Therefore their 
6 becomes 

(5 = 3.704 X 10-2 Z^i^Lo , (30) 

''"304 

and the radiation pressure due to A304 is 

' " 3A304nHe//AC304(l + 0.9,5)0-97 ■ ^ > 

In these models, radiation pressure is dominated by Lyman-a and Hell Lyman-a . 



5. Results 



The mean intensity, at threshold Ji^, of the UV background is calculated at Zon and at 
z = 2.5 for each model. The mean intensity is normalized at z = 2.5 by: model A, 47r = 
Att X 10-21 ergs cmr'^ s'^ Hz'^; model B, = 10 x J^; and model C, = 100 x J^. 
For each model, and Zq^, the maximum cloud size before the cloud becomes unstable to 
radiation pressure {Pr > -fg) is calculated using the photoionization code CLOUDY. Finally 
a model is made of the cloud at z = 2.5 scaling the density and radius of the maximum 
permissible cloud at Zq^ and the neutral column density is recorded. The density at Zon is 
varied and the procedure repeated. The results of our calculations are shown in Figures ?? 
and[6l The area above and to the right of each curve is not allowed because of the instability 
due to radiation pressure at Zon- 
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Expansion Method ^Hl {cm 




5 


10^7 


Isothermal 


10 


10^6 




20 


10^5 




5 


10^6 


Adiabatic 


10 


10^4 




20 


1012 



Table 2: Results for isothermal and adiabatic expansion of Lyman-a clouds. 

There are two distinctive features in each model: a linear or nearly linear phase where the 
maximum permissible cloud size does not vary appreciably with increasing neutral hydrogen 
column density; and a vertical phase where the cloud size changes dramatically and the 
neutral hydrogen column density remains nearly constant despite changing the density of 
the cloud. It is noted that the appearance and location of the vertical phase is independent 
of the normalized value of the mean intensity, however the location in neutral hydrogen 
column density of the vertical phase is dependent upon the redshift chosen for Zon and the 
scaling law for cloud density. This is shown in Table [2j 

Metal free Lyman-a systems are seen with neutral hydrogen column densities approaching 
IQI^ cm~2 and are not seen with Nhi > lO^^'^ cm~'^. We expect that the Lyman-a clouds 
expansion is not strictly isothermal nor adiabatic, but lies in between these two extremes. 
Therefore the instability because of radiation pressure indicates that Zon < 10 and the lack 
of metal free clouds with Nhi ^ lO^^'^ cm~'^ may indicate that 5 > > 10. 

With more observations of Lyman-a forest clouds, definitive upper limits on the size of 
the clouds and on neutral hydrogen column density would decrease the uncertainty in the 
magnitude of the background UV flux and the epoch of Zq^ for QSOs. 

The author gratefully acknowledges J. A. Baldwin, E.R. Capriotti, G.J. Ferland, J.V. Vil- 
lumsen and R.J. Weymann for helpful dscussions and comments on the manuscript. This 
work was supported in part by National Science Foundation Grant AST 000-00. 
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Redshift History of Ionizing Flux 
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Fig. 4. — The ratio of UV flux at redshift z to UV flux at z = 2.5 (at threshold) vs. redshift 
to redshift z = 4, showing that J^j, is nearly constant for all three models in this redshift 
range. 



28 



27 



26 



E 25 



24 



23 



22 



21 



Maximum Stable Diameter, Isothermal Expansion 





1 


1 


1 


1 


1 


1 


1 












O O 




> o 














o 


o 










O 


rP 


o 


o 








□ 






o 















□ 


o 


o 












□ 


o 


















< o " 












H □ 


o 




















^0" = IS 


O 














Zon = 20 


□ 

















11 


12 


13 


14 


15 


16 


17 1 



log N (cm^) 



Fig. 5. — Maximum stable diameter (in cm) vs neutral hydrogen column density (in cm ^) 
for isothermally expanding clouds. 47rJ^j,(2.5) = An x 10~^^ ergs cm~^ s~^ Hz~^. 
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Fig. 6. — Maximum stable diameter (in cm) vs neutral hydrogen column density (in cm ^) 
for isothermally expanding clouds. 47r 7,^^ (2.5) = 47r x 10"^'' ergs cm~^ Hz~^. 



